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Abstract. - We introduce a model for innovation-, evolution- and opinion dynamics whose 
spreading is dictated by unanimity rules, i.e. a node will change its (binary) state only if 
all of its neighbours have the same corresponding state. It is shown that a transition takes 
place depending on the initial condition of the problem. In particular, a critical number of 
initially activated nodes is needed so that the whole system gets activated in the long-time 
limit. The influence of the degree distribution of the nodes is naturally taken into account. For 
simple network topologies we solve the model analytically, the cases of random, small-world 
and scale-free are studied in detail. 



Introduction. In general, the discovery or emergence of something depends on the 
combination of several parameters, all of them having to be simultaneously met. One may 
think of economy, where the production of a good depends on the production or existence of 
other goods (e.g. to produce a car one needs the wheel, the motor and some fioritures). In 
return, this new discovery opens new possibilities and needs that will lead to the production 
of yet new goods (e.g. the simultaneous existence of the car and of alcohol directly leads to 
the invention of the air bag). This auto-catalytic process is a very general process [1-4] and 
obviously applies to many situations not only related to innovation, but also to evolution, 
opinion formation, food chains etc. One may even think of the dynamics of scientific ideas, 
music genres, or any other field where the emergence of a new element possibly leads to new 
combinations and new elements. This feedback is responsible for the potential explosion of 
the number of items, such as observed e.g. in the Cambrian explosion). This "explosion" has 
been shown to be identical of a phase transition in a Van der Waals gas [1] . After mapping the 
above catalytic reactions onto a network structure, where nodes represent items and directed 
links show which items are necessary for the production of others, it is tempting to introduce 
a unanimity rule (UR) : a node on the network is activated only if all the nodes arriving to it 
through a link are activated. Surprisingly, the dynamics of such an unanimity rule, that is a 
straightforward generalization of the usual majority rules of opinion dynamics [5-9], is poorly 
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Fig. 1 - First two steps of UR starting from an initial network of 7 nodes, 2 of them being activated. 
Initially there is only one node among the non-activated nodes that satisfies the unanimity rule. It 
gets therefore activated at the first time step. At that time, there is a new node whose 2 incoming 
links come from activated nodes. It gets activated at the second time step. It is straightforward to 
show that this system gets fully activated at the fourth time step. 

known [1]. In this Letter, we are interested on the spreading dynamics of the UR and try to 
understand which conditions have to be satisfied so that the network gets fully activated in 
the long time limit. 

Unanimity model. - The network is composed of N nodes related through directed links. 
Each node exists in one of two states: activated or inactivated. The number of nodes with 
indegree i (the indegree of a node is defined to be the number of links arriving to it) is denoted 
by Ni and depends on the underlying network structure. Initially (at t = 0) there are Aq nodes 
which are activated, among which A i have an indegree i. In general, the number of nodes 
of type i, activated at time t, is A it . It is also useful to introduce the quantities = and 
a>i,t — which are the proportions of nodes of type i in the network (indegree distribution) 
and the probability that such a node i is activated, respectively. a t = 4f = J2i n i a i,t is the 
fraction of activated nodes in the whole network at time t. The unanimity rule is defined as 
follows (see Fig.l). At each time step, each node is considered. If all the links arriving to a 
specific unactivated node i originate at nodes which are activated at t — 1, i gets activated 
at t. Otherwise, it remains unactivated. The process is applied iteratively until the system 
reaches a stationary state, characterized by an asymptotic value a^. In the following, we are 
interested in the relation between and arj, i.e. what is the final occupation of the network 
as a function of its initial occupation on a specific network. Let us mention the fact that 
each node may be produced by only one combination of (potentially many, depending on the 
indegree) nodes. This is a modification of the model of Hand ct al. [1], where more than one 
pairs of (two) nodes could produce new elements and will lead to a different equation for the 
activation evolution, as shown below. The dynamics studied here implies that nodes with a 
higher indegree will be activated with a probability smaller than those with a smaller indegree 
(because the former have more conditions to be fulfilled). 

Master equation. Let us now derive an evolution equation for Ai tt and A t . To do 
so it is helpful to consider the first time step and than to iterate. There are initially Aq 
activated nodes, = A Ni/N of them being of indegree i on average (the activated nodes 
are randomly chosen in the beginning). The ensemble of A^o nodes is called the initial set of 
indegree i. By construction, the probability that i randomly chosen nodes are activated, is 
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dp (i is an exponent). Consequently, the average number of nodes with indegree i and who 
respect the unanimity rule is Niclq while the number of such nodes that are not yet occupied 
is 

A i>0 = (Ni - iViaoKo, (1) 

and, on average, the total number of occupied nodes with indegree i evolves as: 

A itl = A i<0 + A it0 . (2) 

Let us stress that we have implicitely assumed that there are no indegree correlations between 
neighboring nodes in order to derive EqfT] At the next time step, the average number of 
nodes with indegree i, who respect the unanimity rule and who are outside the initial set is 
(Ni — Nida)a\. Among those nodes, Ao^ have already been activated during the first time 
step, so that the average number of nodes who get activated at the second time step is: 

A i ,i = (JV i -JV i ao)(ai-4). (3) 

Note that Eq[3]is valid because no node in Aj,i also belongs to A^q. This is due to the fact 
that each node can only be activated by one combination of i nodes in our model, so that no 
redundancy is possible between A^i and A^q. By proceeding similarly it is straightforward 
to show that the contributions A^t read 

Ai, t = (Ni - N iao )(ai - 4^), (4) 

with a_i = by convention. The number of activated nodes evolve as 

A iit+1 = M,t + A iit . (5) 

By dividing by Ni, one gets a set of equations for the proportion of nodes € [0, 1]: 

a iit+ i = a i<t + (1 - a )(a\ - a^), (6) 

where the coupling between the different proportions a^t occurs through the average value 
at — J2i n i a i,ti as defined above. Finally, by multiplying by the indegree distribution rii and 
summing over all values of i, one gets a closed equation for the average proportion of activated 
nodes in the network that reads 

q t+ i = a t + (1 - ap) rii(a\ - a l t _i). (7) 

i 

Let us stress that EqJ7]is non-linear as soon as Ni ^ 0, i > 1. Moreover, it is characterized 
by the non-trivial presence of the initial condition ao in the right hand non-linear term and 
is therefore highly non-local in time. Eq[7] explicitly shows how the indegree distribution ri^, 
affects the propagation of activated nodes in the system. 

Theoretical results. - In this section, we focus on simple choices of n.^ in order to appre- 
hend analytically the behavior of Eq(7] The simplest case is rii = 1 for which Eq[7] reads 

Ot+i = a t + (1 - a )(a t - a t -i). (8) 
This equation is solved by recurrence: 



ai = a + (1 - a )a 
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Fig. 2 - Numerical integration of Eq |12l In the left figure, we plot the time evolution of at, together 
with the theoretical prediction 1141 In the right figure, we plot the relation Ooo(ao) that obviously 
shows a transition at ao = 1/2. Eq |14l is in perfect agreement with the numerical integration of Eq|12l 
(line). The agreement with the numerical simulations of the model itself is also excellent. 



and in general 



a-2 = a + (1 - a )a Q + (1 - a )[a + (1 - a )a - a ] 

= oto + (1 - a )a + (1 - a ) 2 o (9) 



t 

a t = ^(l- ao ) u ao = l-(l-ao) t+1 . (10) 



u=0 

This last expression is easily verified: 
f 



= ^(l-ao)"ao + (l-a )(^( 1 -ao) u a -5](l-an)"ao) 

U— If.— U—0 

t £+1 

= ^(l-a rao + (l- ao ) m ao = ^(l-ao) u ao. (11) 

u=0 u=0 

The above solution implies that any initial condition converges toward the asymptotic state 
0,00 = 1, i.e. whatever the initial condition, the system is fully activated in the long time limit. 
The relaxation to = 1 is exponentially fast ~ e ln ( 1 " a o)t > 

Let us now focus on the more challenging case rii — 1 where all the nodes have an indegree 
of 2 by construction. In that case, Eq. [7] reads 

Ot+i = Of + (1 - Oo)(Of - «t_i)- (12) 

The non-linear term does not allow to find a simple recurrence expression as above. Though, 
a numerical integration of EqH21 (by using Mathematica for instance) shows that the leading 
terms in the Taylor expansion of at behave like 

f+i 

at = ^oj + 0(t + 2) (13) 



thus suggesting that the asymptotic solution is 
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Fig. 3 - aoo(oo) obtained for random networks with L — 2N (left figure). The measured indegree 
distribution is no = 0.f4, n\ = 0.27, ni — 0.27... The No Zero and Zero versions (see text)are shown. 
In the right figure, aoo(fflo) for small-world networks L = AiV short-cuts. When A = 1, the measured 
indegree distribution is no = 0.0, n\ = 0.37, 112 = 0.37, etc. When A = 2, it is no = 0.0, m = 0.14, 
ri2 = 0.27... In both figures, the total number of nodes is N — 10000 and the results are averaged over 
100 realizations of the process. The solid lines are the corresponding numerical solutions of Eq |18l 
evaluated with the empirical values rn. 



This solution should satisfy the normalization constraint < 1, so that it can hold only for 
initial conditions ao > 1/2. This argument suggest that a transition takes place at a c = 1/2, 
such that only a fraction of the whole system gets activated when ao < a c while the whole 
system activates above this value (see Fig. 2). We verify the approximate solution EafTH by 
looking for a solution of the form at = l a ° a (1 + e t ). By insterting this expression into Eq[T2l 
one gets the recurrence relations: 

£f+i = e t + a (l + e t ) 2 - a (l + e*-i) 2 

f-i+i = e t +2a (e t - e t _i), (15) 

where the second line is obtained by keeping only first order corrections in e. In the continuous 
time limit, keeping terms until the second time derivative, one obtains 

(1 - 2a )d t e t + 1/2(1 + 2a )dfe t - 0, (16) 

whose exponential solutions read e t = e~ xt with 

This is a relaxation to the stationary state only when ao < 1/2, thereby confirming a 
qualitative change at a c = 1/2. 

Some network topologies. - Let us now focus on more reasonable topologies and compare 
the results obtained from Eq. [7] with numerical simulations of the UR. We focus on three types 
of networks, purely random networks [10], small- world like networks [11] and Barabasi- Albert 
networks [12] (growing networks with preferential attachment). The excellent agreement with 
Eq. [7] suggest that the formalism should apply to more general situations as well. The 
random network was obtained by randomly assigning L directed links over N nodes. The 
small-world network was obtained by starting from a directed ring configuration and than 




Fig. 4 - (left figure) distribution of indegrees for the scale-free network, rrii = m+i (this shift being 
made to plot no in log-log scale). The solid line is the theoretical asymptotical prediction k~' 1+1 ' p ', 
for p — 1/2. The total number of nodes is N = 10000, results are averaged over 100 realizations of 
the process, (right figure) aoo(ao) for the same network. In the No Zero version, one observes that 
the all the nodes are finally activated whatever the initial condition. The solid lines correspond to 
numerical values of Eg 1181 for the empirical values n;. 



randomly assigning L directed links (short-cuts) over the nodes, i.e. the total number of links 
in that case is L + N (The network drawn in Fig.l is such network with N = 7 nodes and 
L = 3 short-cuts). Let us note that the small- world network can be viewed as a food chain 
with a well-defined hierarchy between species together with some random short-cuts. In that 
case, UR can be interpreted as an extinction model (if all the species that one species eats 
go extinct, this species will also go extinct). The Barabasi- Albert network was built starting 
from one seed node and adding nodes one at a time until the system is composed of N nodes. 
At each step, the node first connects to a randomly chosen node and, with probability p, it 
re-directs its link to the father of selected node. This method is well-known to be equivalent 
to preferential attachment and to lead to the formation of fat tail degree distributions k u : 
with v= 1 + 1/p, [13]. 

Once the underlying network is built, we randomly assign a N active nodes to the network 
and apply the unanimity rule. The evolution stops once a stationary state is reached. The 
asymptotic value doo is averaged over several realizations of the process (on several realizations 
of the underlying network). In the small- world network, each node receives at least one 
incoming link. This is not the case for the random- or the BA networks, for which one has to 
discuss the ambiguous dynamics of nodes with zero incoming links. Two choices are possible. 
Either these nodes can not be activated in the course of time, because they are not reached 
by any other node (No Zero version), or all of them are get activated at the first time step, 
thereby assuming that their activation does not require any first knowledge (Zero version). The 
choise is a question of interpretation. The two versions are associated to different evolution 
equations: 

oo 

at+x = at + (1 - ap) nj{a\ 

z=l 
oo 

<H+i = at + (1 - a )22 n i( a t 

i=0 

and leads to quite different behaviors (Figs. 3 and 4). In the case of small- world networks, the 



cit-i) N° Zero version 

aj_i) Zero version, (18) 
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above equations are obviously equivalent. To compare the simulation results with Eq |18l we 
also measure the indegree distributions of the networks rii generated during the simulations 
and integrate Eg J18I with these empirical values. The agreement is excellent, except close 
to the transition points where finite size effects are expected. It is worth noting that the 
importance of nodes with zero incoming links is much higher in (growing) Barabasi- Albert 
like networks (Fig. 4), so that the difference between the two versions is quite pronounced, as 
expected. Let us also mention that Eq.14 has been successfully verified for a random network 
where the indegree of each node is exactly 2 (712 = 1), as shown in Fig. 2. 

Discussion. - In this Letter, we have introduced a simple model for innovation whose 
dynamics is based on Unanimity Rule. It is shown that the discovery of new items on the 
underlying network opens perspectives for the discovery of new items. This feedback effect 
may lead to complex spreading properties, embodied by the existence of a critical size for the 
initial activation, that is necessary for the complete activation of the network in the long time 
limit. The problem has been analyzed empirically on a large variety of network structures and 
has been successfully described by recurrence relations for the average activation. Let us stress 
that these recurrence relations have a quite atypical form, due to their explicit dependence 
on the initial conditions. Moreover, their non-linearity makes them a hard problem to solve 
in general. Finally, let us insist on the fact that Unanmity Rule is a general mechanism 
that should apply to numerous situations related to innovation, opinion dynamics or even 
species/population dynamics. To be consistent with our own work, we also hope that this 
paper will trigger the reader's curiosity and, possibly, open new perspectives or research 
directions... 
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